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SUMMARY 


This is a working report prepared in order to document early 
results of research on the stability of laminar boundary 
layers. The report shows that constitutive equations for 
a structured continua may be derived by the simple technique 
of reinterpreting velocity in the conventional stress to 
rate-of- strain relationship so as to account for effects of 
particle rotation. 

The report also demonstrates that accounting for particle 
structure even at a molecular level makes the fluid visco- 
elastic with the ability to propagate vector waves. Finally, 
it is shown that particle structure modifies the basic 
stability equation for the system, which in turn would alter 
values for critical Reynolds number. 

It therefore has been demonstrated that conventional fluids 
such as air actually are visco-elastic, and that the Navier- 
stokes equations do not furnish a proper mathematical model 
for the investigation of hydrodynamic stability. 
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I. INTRODUCTION 

The location of boundary- layer transition along a surface 
is one of the factors necessary for the prediction of aero- 
dynamic performance in highspeed flight. This is apparent 
when one considers that the condition of the boundary layer 
(i.e., whether laminar or turbulent) determines heat 
transfer, skin- friction drag, base-pressure drag, and 
effectiveness of control surfaces that are partially immersed 
in the body boundary layer. Furthermore, effects of the 
boundary layer on various types of "pressure" drags and 
control-surface effectiveness become more important as the 
flight speed increases, due to the fact that boundary- layer 
thickness increases with Mach number (at constant Reynolds 
number) . 

Around 1950, there was little quantitative information 
available regarding the effects of various parameters on 
boundary- layer transition, especially in supersonic flow. 
Research on subsonic flow had shown that the transition 
Reynolds number is influenced by surface temperature, sur- 
face roughness, pressure gradient, surface curvature, and 
free-stream turbulence. Investigations in transonic flow 
indicated that shock-wave/boundary- layer interaction also 
influences transition. However, only the qualitative 
effects of these variables were known; little quantitative 
information was available about the effect of these variables, 
either singly or interrelated. Quantitative information 
about boundary layer stability resulted from the theoretical 
work of Tollmien (1936) and Schlicting (1935) . For years 
later, no experimental evidence was found to confirm the 
laminar oscillations predicted by Tollmien-Schlicting theory. 
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Finally, in 1947, the existence of Tollmien-Schlicting 
oscillations was confirmed by the experiments of Schubauer 
and Skramstad. However, this theory and its extensions only 
yield the Reynolds number above which disturbances may 
become sufficiently amplified to cause turbulence. The 
theories are useful, however, for indicating relative 
significance of various parameters on transition. 

Because of the many factors that influence boundary- layer 
transition, the experimental study of the individual effect 
of any one of these variables on transition is difficult. 

For example, the effects of free-stream turbulence in wind- 
tunnel studies tend to mask out the effects of other 
variables. This is especially true of supersonic flow, since 
little is known about turbulence in supersonic flow or its 
effects; furthermore, the turbulence characteristics of most 
of the supersonic tunnels have not even been completely 
determined to this day. 

The theory of stability of laminar flows was first formulated 
mathematically by Rayleigh (1887). Successful solutions for 
the laminar boundary layer on a flat plate were obtained by 
Tollmien and Schlicting, as already noted. The results of 
this theory indicate that instability in the laminar boundary- 
layer originates in small disturbances that are either damped 
or amplified depending on their frequency and the Reynolds 
number. This selective amplification is predicted only if 
the Reynolds number is greater than a certain value, called 
the "minimum critical Reynolds number". This Reynolds 
number is always found to be somewhat smaller than the 
transition Reynolds number, partially due to the fact that 
regular disturbances in the laminar boundary layer must have 



time to be amplified to a magnitude sufficient to somehow 
produce the irregular oscillations characteristic of 
turbulence . 

A somewhat different theoretical approach was suggested 
by Taylor (1936). The theory of Taylor, originally 
thought to be in opposition to that of Tollmien and 
Schlicting, involves the supposition that the free stream 
outside the boundary layer always carries at least a slight 
turbulence, and this turbulence imposes local adverse 
pressure gradients on the laminar boundary layer. These 
local adverse gradients are thought to cause local boundary 
layer separation (insipient separation), with subsequent 
transition. 

The apparent difference between these two theories was 
resolved by the experimental studies of Schubauer and 
Skramstad. In connection with a study of wind-tunnel tur- 
bulence, they observed the regular laminar boundary- layer 
oscillations predicted by Tollmien and Schlicting in a 
wind tunnel having an exceptionally low turbulence level. 
The predicted ranges of damping and amplification were 
closely confirmed. However, it was found that the oscilla- 
tions and their damping or amplification to eventual tur- 
bulence could be detected only when the magnitude of the 
free-stream turbulent fluctuations is less than about 0.001 
of the mean velocity. For larger magnitudes of the free- 
stream turbulence, the boundary- layer oscillations are 
difficult to identify because of the near coincidence of 
their appearance and the point of transition to turbulence. 
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With the aid of the Schubauer-Skramstad experiments, the 
Tollmien- Schlicting theory and the Taylor theory can be 
reconciled, and a unified picture of the mechanism leading 
to transition in the boundary- layer can be obtained. In 
the case of low free-stream turbulence < 0.001 j , 

any small disturbances that are initially present may be 
due to either internal or external disturbances such as 
surface roughness, vibration, noise, or free-stream tur- 
bulence. The amplified oscillations themselves cannot be 
classified as turbulence because of their regularity; it 
is only when these oscillations or waves become large enough 
to roll up into random eddies that turbulence is initiated. 
The eddies, formed close to the solid surface, break away 
and dissipate in an irregular turbulent motion. The process 
is well illustrated in the interferometer photographs of 
Eckert (1951) , and is somewhat analogous to the growth and 

breaking of waves on an air-water interface. In the case 

/u * \ 

of a high free-stream turbulence level I > 0.001) , the 
free stream fluctuations are large enough to produce local 
adverse pressure gradients that cause local separation with 
subsequent transition in the manner described by Taylor's 
theory . 

In recent years, numerous revisions and extensions to the 
Tollmien-Schlicting theory have appeared. Lin (1946) has 
made some revisions and mathematical clarifications. 
Schlicting and Ulrich (1942) and Hahneman, Freeman, and 
Finston (1948) extended the theory to flows with pressure 
gradient and found that a positive pressure gradient 
destabilized while a negative pressure gradient stabilized 
the boundary layer. 
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Although the theories of the Tollmien-Schlicting type permit 
prediction of the effect of various factors on laminar 
boundary- layer stability and prediction of the Reynolds 
number below which disturbances are not amplified, they give 
no indication of the transition Reynolds number. As already 
noted, presumably the original small distrubances must have 
time to be amplified to the point where they become unstable. 
While the theory permits calculation of the initial amplifica- 
tion rate, it is strictly applicable only when the disturbances 
are infinitesimally small. That this amplification takes 
considerable time and distance is exemplified by a comparison 
of the minimum critical Reynolds number with the transition 
Reynolds number. For incompressible flow along a flat plate, 
the minimum critical Reynolds number had been calculated 
as 60,000; this value was confirmed by Schubauer and 
Skramstad, who also found the Reynolds number range of the 
transition region to be 2,800,000 to 3,900,000. For this 
case, then, there is about a sixty-fold difference between 
the Reynolds numbers at which amplification begins and at 
which turbulence begins. 

Liepmann (1945) has made an interesting attempt to compute 
the transition Reynolds number using the maximum initial 
rates of amplification given by the Tollmien-Schlicting 
theory. He estimates that the transition to turbulence 
will begin when the apparent shear stress of the amplified 
disturbances becomes equal to existing laminar shear stress. 

The resulting relationship gives the transition Reynolds 
number as a function of the magnitude of the initial disturb- 
ance having the frequency of maximum amplification. However, 
this latter quantity is dependent on the magnitude and 


<1 
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frequency of the free-stream turbulence, of the surface 
roughness, and/or of noise and vibration; the quantitative 
dependence of this factor on these variables is not known. 
Therefore, Liepmann's calculation can’t be compared directly 
with experiment; however, for reasonable values of the initial 
disturbances, the predicted ratio of transition Reynolds 
number to minimum critical Reynolds number is of the same 
order as that found experimentally. 

Emmons (1951) has evolved a theory of transition based on 
the idea that every point on a body is turbulent part of 
the time, and the fraction of the time that a given spot 
is turbulent is derived from probability considerations. 

The results of this theory, while not enabling prediction 
of the absolute value of transition range of Reynolds 
numbers, give the variation of skin friction or heat- transfer 
rate in the transition region and the effect of flow geometry 
(flat plate, cone, swept wing) on the transition. 

Basically, all work in recent years on boundary layer 
stability has followed the Tollmien-Schlicting pattern, 
based on better and better analysis of the Navier- Stokes 
equations as descriptive of the system. More recently, 

Kistler (1969, 1971) has reviewed existing literature and 
assembled clues to the effect that the Navier-Stokes equa- 
tions themselves may not provide an adequate specification 
of the stability of the real fluid system being modeled. 

Based on a review of literature with a wider scope than 
conventional boundary- layer stability theory, a coupled- 
dynamics model (Kistler, 1971) based on the work of Dahler 
(1959) and associates has been suggested as more representa- 
tive of the flow of diatomic gases. This report documents 
an initial investigation of the stability of that model. 
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The stage is set by noting it has been recognized for many 
years that internal degrees of freedom of molecules play 
important roles in determining microscopic thermodynamical 
properties of matter. However, the classical theories of 
fluid mechanics (i.e. Navier-Stokes equations) do not include 
explicitly the kinematical effects of these same internal 
degrees of freedom, apparently because it is not widely 
recognized that in some situations these kinematical effects 
may play an important role in stability of the flow and in 
nonequilibrium processes involving momentum transport. 

This report presents first a simple derivation and interpre- 
tation of the constitutive equations for fluid species with 
a finite minimum scale and the ability to transfer intrinsic 
angular momentum between molecules. Diatomic or rough-sphere 
monatomic fluids would fall in this category. 

In work on the theory of structured continua, Dahler (1959), 
Dahler and Scriven (1963) , Condiff and Dahler (1964) , and 
others derived these same equations by a more elegant but 
involved procedure, and showed that internal angular 
momentum is associated with configurational and kinetic 
structure of the continuum "particle". Condiff and Dahler 
assign a continuous spin field to the rotation or spin of 
molecular subunits. The interaction of this internal spin 
with fluid flow is described by antisymmetric stress, while 
couple stress accounts for viscous transport of internal 
angular momentum. Then, with constitutive equations appro- 
priate to a linear, isotropic fluid they obtain generalized 
Navier-Stokes equations for the velocity and spin fields. 
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The macroscopic spin field is defined by the local average 
rotational velocity of the molecules, and is a kinematical 
representation of internal angular momentum of spinning 
molecules. This spin field is dynamically coupled to the 
fluid velocity by means of the collisional interactions of 
the translating and rotating molecules. From a kinetic 
theory point of view, this coupling is a result of noncentral 
intermolecular forces. Condiff and Dahler note that Born 
was evidently the first to suggest that antisymmetric stress 
is the result if this internal spin is not "synchronized" 
with the vorticity. 

This document utilizes this background which establishes the 
coupled-dynamics math model, and goes on to show that since 
including the effects of molecular structure has coupled the 
translational and rotational equations of motion, the system 
now can actually support the propagation of vector waves 
(transverse waves) within the interior of the fluid. These 
waves are waves of angular momentum (or angular velocity) , 
and will be called "spin waves". 

In classical theory of hydrodynamic stability of viscous 
fluids such vector waves cannot exist. Specifically their 
existence has been precluded by the neglect of molecular 
substructure . 

Finally, in this report it is shown that under assumptions 
comparable to those leading to the classical Orr-Sommerfeld 
equation, the analogous equation accounting for substructure 
contains an exponential term which could lead to instability 
of translational perturbations. 
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II. Constitutive Equations 

Conservation equations for a structured continua have been 
developed and presented in a number of journal articles and 
texts. In the notation used by Condiff and Dahler (1964), 
the equations for conservation of linear momentum and intrinsic 
angular momentum of a steady incompressible flow are 

+ u*Vu^ = pf + V*S II - 1 

l * ~"1 ~ aa A 

— + u* V£ = pG + V*C + A . 1 1 - 2 

Constitutive equations for the stress and couple stress 
tensors were developed by Dahler (1959) , and Dahler and 
Scriven (1963). It is found here that those same constitu- 
tive equations can be obtained by a much easier procedure 
if the relative velocity between two particles in the fluid 
is properly interpreted. 

From the theory of structured continua, the intrinsic angular 

A 

momentum $L is taken as 



l = I*a) . 1 1- 3 

o 

N\ AA /s 

where the tensor field I = I(r,t) simulates the average 
moment of inertia of the molecules. However, if the molecules 

M 

are rigid then I can be regarded as a constant dyadic, and 
if the fluid is isotropic* then 

AA 

I = I U 1 1 - 4 


& 

This is a critical assumption common to all models todate. 
It allows relating an average angular velocity to voritcity 
as opposed to an angular momentum relationship. 
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Here U is the unit dyadic and I is a constant scalar. 

A'N 

See Appendix A for second or alternate interpretation of I . 


If consideration of the constitutive relations for . stress is 
restricted to the linear Stokesian regime, the stress-strain 
relationship is found by noting that the velocity field which 
determines the velocity gradient is not strictly defined by 

A 

the translational velocity field u , but instead it is 

A 

represented by a local velocity U which may be slightly 

A 

greater or less than u due to rotation of the molecule. 


If the molecules are assumed rigid, then the tangential 
velocity of an atom due to rotation of the molecule is 

A A /N 

Rxto , where R is the average radius of a molecule in the 
flow. The atom then moves with an absolute velocity U 
which is 


U 


u + Rxco 


o 


II- 5 


and 


VU = Vu + V(Rxto ) 

o 


Vu + [VRxoj^] - [V(jJ q xR] . 


1 1 - 6 


However , 


VR = U • 


II- 7 


Thus , 


VU = Vu + (Uxu> ) + [RxVoj 1 . 

o L o 


II- 8 
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Even though the spin field may vary from particle to 

particle as does the translational velocity u , the spin 

A 

vector a) is constant for any one particle if the particle 
is rigid. Hence Vw q = 0 for an individual particle and 
does not influence the absolute velocity at the edge of the 
particle. Therefore, the absolute velocity gradient is just 

VU = Vu + fljxu> ) . 1 1 - 9 


The last term clearly is skew-symmetric, and may be combined 
with the skew part of the gradient of translational velocity. 
Thus 

VU = (Vu) (s) + (Vu) (a) + (l)xto o ) 

= j (Vu + uV) + j (Vu - uV) + (U*w o ) 11-10 

= j (vG + Gv) - j U x<i (vG - Gv)^> + (Uxu o ) , 

where <^^-(Vu-uV)^> is the pseudovector 
translational velocity dyadic. 

Now the pseudovectors of Vu and uV a 

<(Vu) = Vxu 

and 

<uV) = -Vxu 


of the skew part of 


re 


11-11 
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hence 


< 


1 

7 



V*£i . 


11-12 


Thus , 

VU = i (Vu + uV) - j U x (VxG) + (U x w o ) 

= j (Vu + uV) - 2 U x [VxG - 2w q ] 11-13 

= i (Vu + uV) - i dual w => V + co 

l l ~ mn ~ mn 

c 

where w is a relative angular velocity due to the intrinsic 

~ A 

spin field uj 0 not being "synchronized" with the vorticity 
Vxu. it will be seen later that oj is like a relaxation 
parameter. 


In the simple case where it is assumed that the relationship 
between stress and rate of strain is linear and isotropic, 
the constitutive equation is written as 


T. . = D. . [V + a) ] 

i j x j mn mn ~mn J 


11-14 


Here V mn are the components of the symmetric velocity 
gradient dyadic, are the components of the skew part, 

and 


D.. = a6 . . 6 +06. 5. +y<$.6. , 11-15 

i]mn X] mn xm ]n in ]m 


pq 


where 6 


is the Kronecker delta. 
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The D. . are components of the viscosity tensor, and this 
tensor is asymmetric, so it too may be written in terms of 
symmetric and skew parts. Let a = A ; 0 = y + 6 ; and 

Y = y - 3 . See discussion in Appendix B. 


Then , 


T. . = r A6 . .6 + y (6 . 6 . + 6.6.) 

13 |_ 13 mn 1m ]n in 3m 

+ 6(6.6. - 6 . 6 . )1 (V +03 ) 

1m ]n m 3m J mn ~mn' 


11-16 


Performing the indicated multiplication gives 


id 


A6 . .6 V + y (6 . 6 . + 6 . 6 . )V 

13 mn mn v 1m 3n in 3m mn 


+ 6 (6 • 6 • - 6 . 6 . )V + A6 . .6 o) 

1m jn in 3m mn 13 mn~mn 


+ y(6. 6. + 6 . 6 . ) o) +6(6. 6. - 6 . 6 . ) 03 

1m jn in jm ~mn , 1m jn m jm ~ 


mn 


11-17 


Since the product of a symmetric and a skew- symmetric tensor 
is zero, the viscous stress tensor , reduces to 


T. . = A6..6 V +y(6. 6. + 6. 6.)V 

13 ij mn mn 1m 311 in 3m mn 


+ A6..6 03 +6(6. 6. - 6 . 6 . ) 03 

13 mn~mn 1m jn in 3m ~mn 


A6 . .V + 2yV. . + A6 . . 6 03 + 6e, . . £, 03 

ij mm ij 13 mn-mn kij kmn~mn 


11-18 
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But 6 mn is symmetrical, so 


6 03 = 0 . 

mn~mn 


11-19 


Hence, the viscous stress tensor is just 

T = xUv*u + y [Vu + uV] + B dual u> , 11-20 

and the full stress tensor is 

/\A , , .a A A A A A 

S = -pU + AUV*u + y [Vu + uV] + 8 dual(Vxu - 2a> ) . 

11-21 


The coefficient y is the usual coefficient of shear 
viscosity, and A is second coefficient of viscosity. The 
third coefficient 3 is one half one., what Condiff and 
Dahler have recently called a coefficient of vortex viscosity. 


Contracting Eqn. 11-18, identifying p with the average 
hydrostatic pressure*, and identifying k = 3A. + 2y as the 
usual coefficient of bulk viscosity** allows the stress 
tensor to be written as 

's = U(~P + j kV»u) + y [(Vu + uV) - y Uv.u] 

+ 6 dual [Vxu - 2a3 Q ] 11-22 


Note that this is not really valid for problems involving 
finite disturbances, high frequencies, and relaxation 
effects. The analysis should be generalized further. 

** 


See Appendix B 
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This is the same result derived by Dahler by a more elegant 
but involved procedure employing kinetic theory, whereas 
here it is obtained by a fairly simple procedure linked 
closely with intuition and without the need to resort to 
kinetic theory. 

Condiff and Dahler take the couple stress tensor to be 
symmetric and to depend only upon the symmetrized spin 
gradient tensor. Hence, by a similar procedure to that 
above, the couple stress tensor is 

'c = v U V’w + v„[Vcj + a) V - \ Uv*w ] 11-23 

1 o 2 L o o 3 o 

It is clear from Eqn. 11-22 that the skew part of the stress 
tensor is 


/ S (a) = 3 dual (Vxu - 2u ) , 

O 


and the pseudovector of this tensor is 


A = 23 (Vxu - 2 oj q ) , 


11-24 


11-25 


which appears in Equation 1 1 - 2 . 
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III. Spin Waves 

Equations 11-22 and 11-23 are the constitutive equations for 
a structured continua, and represent the relationship between 
stress and rate-of-strain when the rigid-body rotational 
degree of freedom is accounted for. When these expressions 
are substituted into the kinematical equations for balance 
of linear momentum and balance of intrinsic angular momentum, 
i.e., Equations I I - 1 and II-2, then the resulting equations 
for the fluid system are* 

p|^jr + u*Vuj = -Vp + |k + j p - gjw*u 

+ (y + g)V 2 u + 2 3 V x co q 1 1 1 - 1 

pi 

+ 2g ( Vxu - 2oo q ) 1 1 1 - 2 


3 co 

O /\ 

tytt— + u*Vw„ 
3 1 0 


= ( V 2 + I V l) 


VV*cj 0 + ( v 2 + v i )V‘ 


ur 


This is a set of nonlinear partial differential equations 
representing the dynamics of the system, and demonstrates 
the coupling between the translational field u and the 
molecular spin field w 0 . The vortex viscosity serves as a 
coupling coefficient because if g = 0 then the spin field 
does not influence the translational motion. 


See Kistler (1971) for references on the derivation of 
these equations. 
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First note that for the uncoupled case (0 = 0) the equations 
are 

p ^jpjr + u*Vuj = -Vp + |k + yjvv*u + pV 2 u III-3 

(v 2 + i v^VV-^ ♦ (v 2 ♦ vpv 2 £ 0 . 

1 1 1 - 4 


' 3 % — 

jt~ + u ' Va) o. 


These are of the form 



1 1 1 - 5 


1 1 1 -6 


which are parabolic and representative of diffusion phenomena. 
Arguments applied to Equation III -4 produce the classical 
result of conservation of vorticity, and time rate of change 
of vorticity via diffusion is obtained from the curl of 
Equation 1 1 1 - 3 . By Helmholtz's theorem, 

u = curl A + grad B , 1 1 1 - 7 

where A is some vector and B is a scalar, it can be shown 
from Equation 1 1 1 - 3 that there exists a scalar wave equation 
for the system, representing wave propagation of potential B . 
However, the equation resulting from the vector potential 

A 

A is again a diffusion equation, so that the system will not 
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support propagation of transverse (i.e. vector) waves. 

The scalar wave equation from this procedure is just the 
equation for a longitudinal wave, and describes propagation 
of sound in the fluid. See Morse and Feshbach (1953). 

At first glance the presence of coupling terms 28Vxw 0 
and 28 (Vxu- 2to 0 ) does not appear to change this situation. 
Since 

2u> = V*u - 2co 

o 

= t - 2u> 1 1 1 - 8 

and 


28Vxaj^ = 8VxVxu - 2 8 Vxy 


III -9 


where 


VxVxu = VV*u - V 2 u , III-10 

it appears that Equation 1 1 1 - 1 remains parabolic and the 
nature of the diffusion process is only modified slightly by 
the term 28Vx(o. 

Since w is taken to be a molecular rotational relaxation 
parameter and 8 is the associated spin viscosity coefficient 
as discussed in Section II, the magnitude of this coupling 
term is expected to be very small (see Appendix A) . One 
would then conclude that the dynamics of the system are 
adequately revealed by the uncoupled case expressed by 
Equation 1 1 1 - 3 . This equation of course is the Navier-Stokes 
equation, and is classically taken as the equation of motion 
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for the system. From it one obtains the Orr-Sommerfeld equa- 
tion which is used to evaluate the stability of the system. 

At this point, it seems that the skew-symmetric part of the 
stress tensor has not added anything significant to the 
problem. However, it will be shown now that this is not the 
case. Clearly, the skew-symmetric term in Equation III-l 
did couple this equation to Equation III -2. Thus, if (3^0, 
then another mode of energy transfer has been established. 

If the translation and rotation are oscillatory there is 
a potential resonance between the two modes, and even a very 
small coupling term can exchange significant energy over a 
sufficient number of oscillations. 

This is particularly true when the shear coefficient of 
viscosity is small, and then the system is near neutral sta- 
bility, i.e. near transition. Then the dynamics of the system 
could be tremendously influenced by coupling an additional 
degree of freedom. Experience with lightly damped systems 
having large inertia terms indicates that the dynamics of 
these systems can be profoundly influenced by normally trivial 
terms if the system is near neutral stability. This suggests 
there may be a difficulty in numerical computation of the 
stability of boundary layers near transition due strictly 
to the techniques used to suppress numerical instabilities 
since these may actually mask physical instability charac- 
teristics which are being studied. 

3 

Taking ^ of Equation III -2 gives an equation describing 
spin response of the system to a jerk which provides sudden 


c 
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change in torque on a fluid element (recall that the element 
now has dimensions and inertia) : 

/ i v ^ 2 

( v 2 + T v i) vv '7t £ * tv 2 * V v rr 

' A 

+ 23 (Vxu) - 43 ^ . III-ll 

This appears still to yield no further information, but merely 
represent a diffusion equation for the vector quantity 

by virtue of the first term on the left side and the 

second term on the right side. But consider now the coupling 
term (third term on RHS) and see what it does to the equation. 
For simplicity, take the fluid to be steady and incompressible. 
Then using 

V 2 u = VV-u - VxVxG , I I I -12 




Equation 1 1 1 - 1 becomes 


P 



+ 


u • Vu 


-Vp + 


(‘ 


+ 



VV *u 


(3 + y) Vx£ + 23Vxu> 0 , 

III -13 


Taking the curl of Equation III-13 removes the pressure term 
and gives 

P (Vxu) + pVx ( u • Vu) = -(3 + y)VxVx£ + 23VxVx^ Q 

III-14 
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3 a 

Solving for (Vxu) from Equation III-14 and substituting 
into Equation III-ll yields 


Pi 


3 

Tt 


(tT * 2 ' v£ o) 


y 2 * J"l ) VV -^ + ( 111 * V 2 . - V V STT 

3uS ( 

- 43 + 23|-Vx(u*Vu) 


p + 


y) Vx\7x£ 


23 

P 


VxVxw 


1 1 1 - 1 5 


This equation is rearranged to give 



VxVxw Q = VV*w 0 - V 2 w Q , 


1 1 1 - 16 


Since 


1 1 1 - 1 7 


it is clear that this equation has combined features of 

parabolic and hyperbolic , and the hyperbolic character is 
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representative of wave propagation! The parameter being 
propagated is uj Q , which is a vector quantity! Thus, 
this system can actually support the propagation of vector 
or transverse waves, in contrast to all modern classical 
theory about the mechanics of viscous fluids. 


Note that a suitable combination of the third and eighth 
terms in Equation III-16 represents diffusion propagation 
of the vorticity £ because of Equation III-8, as would be 
expected to exist in accordance with classical theory. 
However, note also that a suitable combination of the second 
and eighth terms represent wave propagation of vorticity. 
Finally, note that a suitable combination of the second and 
sixth terms represents a diffusion of the vector quantity 



Consider now the simplified case where the system is still 
coupled (8^0), but let the couple stresses vanish 
(v = v 2 = 0) • Equation III-16 then reduces to 


VxVxto 


0 


p 2 l\ 3 % 


.48 


) d 0 

8t : 


p\ 3w o /p 2 l\ 9 r ~ n - n 

r) rr-[^)n [u - v “o ] 


jg (3 + y)Vxyx£ - Vx[G*Vu] 


III-18 


Replace the spin vector by 


2w, 


Vxu - 2y = V x u - R , 


III-19 


to give 
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III -20 


This still is a vector wave equation in ft , and in £ . 
Note that if the fluid is "relaxed" and R = 0 , then 
Equation III -20 reduces to 


VxVx£ + 





8 

at 


A A 

[u*V5] 


+ 



[u* Vu] 


0 


III -21 


which appears to say that even when the rotational mode is 
relaxed the system still can support a vector wave, but such 

A 

is not the case. Under the condition where R = 0 , and 
with the assumption that v = v 2 = 0 , it is clear from 
Equation 1 1 1 - 2 that the classical case of conservation of 
vorticity has been retrieved from this more general problem. 
That is, 


s\ 

d? 

ar 



A /V 

u*Vc 


0 


1 1 1 - 2 2 


so Equation 1 1 1 - 2 1 reduces to 

VxVx£ + (9-\ II + |£.jvx[u’Vu] = 0 , III- 23 

which is not a vector wave equation, but merely a vector 
analog of the vorticity diffusion equation. 
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Classical sound waves are scalar and are basically waves of 
compression and rarefaction, but the only component of the 
particle displacement that is propagated as a wave motion 
is that lying in the direction of propagation. For this 
reason they are called longitudinal. 

In contrast, the transverse waves have a vectorial character 
in that it is necessary to specify the orientation of the 
displacement in the yz-plane (i.e. a plane normal to the 
direction of propagation). It is clear now that a fluid 
described by Equations 1 1 1 - 1 and 1 1 1 - 2 can support both 
longitudinal and transverse waves, and that they may be 
coupled. This then is at least one mechanism whereby strictly 
2D disturbances can and do become 3D. Except for the obvious 
additional possibility that no experimenter can ever produce 
a completely 2D disturbance, the mechanism discovered here 
is the first answer ever provided for 2D to 3D transfer in 
the linear regime of the problem. 

Longitudinal waves are added algebraically, as scalar quan- 
tities. Transverse waves are added vectorially, with com- 
ponents in two mutually perpendicular transverse directions 
added separately. Now this vector property gives rise to the 
various phenomena of polarization, and so it is reasonable to 
expect that if sufficiently accurate experiments were per- 
formed on the boundary layer region they would indicate that 
the molecules become mechanically polarized due to flow along 
a boundary. 


Before moving on to the question of stability of the system, 
it should be noted that the classical absence of transverse 
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waves, in a fluid is essentially attributed to its inability 
to support shear stress. In saying that transverse waves 
do not exist in a fluid, this of course is speaking of the 
interior of a mass of fluid and beyond reach of any surface 
effect, because it clearly is a matter of common observation 
that transverse waves can exist on the surface of a liquid. 

There is a clear-cut distinction between the longitudinal 
and transverse waves in the system. There is no possibility 
of regarding the longitudinal wave as merely the first com- 
ponent of some more general vector wave whose second and 

/\ 

third components constitute the transverse wave for . 

The two waves are separate solutions for the dynamics of the 
system, and the two waves travel at different velocities of 
propagation. Furthermore, if the solution starts out longi- 
tudinal it will remain longitudinal and if it starts out 
transverse it will remain transverse i_f the coefficients in 
the equation have everywhere the same value. However, if 
the coefficients can vary in a region of space then wave 
reflections will occur and the longitudinal and transverse 
waves may become intermingled. 

In electromagnetic theory, transverse waves are associated 
with the propagation of angular momentum. Transverse waves 
also exist in theory of elasticity, and again are associated 
with angular momentum. However, angular motion in a solid 
is constrained since a given element is not completely free 
to rotate. However, it may exhibit a rotational oscillation 
(Brillioun, 1964) with spring stiffness proportional to 
angular displacement. Counter rotations of adjacent elements 
(Brillioun; and Kistler, 1969) represent a shear phenomena, 
and so transverse waves in a solid are also called shear waves. 
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It is suspected that this is the same type of phenomena 
exhibited by the system under discussion here, and that in 
essence what has been discovered here is that without anything 
like specification of frequency-dependent viscosity coeffi- 
cients, etc. it has been established here that the inclusion 
of molecular characteristics which allow transmittal of 
angular momentum have in effect shown that such a real fluid 
has the properties known generally as visco-elastic! 

Visco-elasticity normally is assumed to exist, via postulates 
about the viscosity coefficients in the constitutive equations. 
Here, in the development of Sections II and III, no postulates 
about visco-elastic or frequency dependent viscosity coeffi- 
cients have been resorted to. All that has been postulated 
is that the fluid does have a structure, and can support 
angular momentum exchange, which of course is well-established 
even for air. Now, it further is well-established that visco- 
elasticity changes the dynamics of the system, and particularly 
the critical Reynolds number obtained.* 


*See a discussion of this in Kistler (.1971) . 
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IV. Stability Analysis 

The governing equations for the system were given as 


[St ♦ H 


■Vp 


+ y u - gjvV’U + (y + 3 )V 2 u + 23VxoS ( . 


IV-1 


pi 


3w 

3T~ + u * Vu) o 


A A 


( v 2 * T v i) vv- “o * (v 2 * v i Jv2 “o 


+ 23(Vxu - 2oa 0 ) 


IV-2 


To show the stability effects of the coupling, it is suffi 
cient to consider a simplified case. Take the flow to be 
incompressible so that V*u = 0 , and let 


A 1 - 

1 

P 

r * 

i 

J y • 

' 6 ) 

> B 1 


A 2 " 

1 

P 

(y + 

3 ) 


J B 2 

" ?r Cv 2 + V 

A 3 = 

23 

P 




5 B 3 

2 3 A 3 

= PT = T ‘ 

Then 

for 

v i : 

= V 2 = 

= 0 

, the 

governing equations 


IV- 3 


9u x a 

St u-Vu = 


.a 


+ A„V u + AVxod 
p 2 3 0 


IV-4 


and 
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3 to 

^ + u • Voo = B _(Vxu - 2uL) = 2B co . IV-5 

dl U 3 U 3 ~ 


Assume 2D parallel flow with the vorticity and spin vectors 
aligned due to mechanical polarization. The governing 
equations in component form then are 


3u 

3t 


+ 



v 


3u 

3y 


3 v 
St 


+ 



v 


3v 

3y 



IV-6 


IV- 7 


3 co 

St 


3 co 
3 


v 


3(0 

Sy 


2B 3 y , 


IV- 8 


where co and y 

~ j 

co Q and (o . 


are the z-components respectively of 


Now make the assumption that velocities and pressure may be 
represented by the sum of a mean term plus a small 
perturbation : 

u=U+u' ; v = V + v' ; p=P+p' IV-9 


The spin field co also is composed of a mean field plus 
perturbations of the same scale as u', v', and p* , as 
well as a relaxation perturbation oo which could be of a 
greatly smaller scale. This would suggest that there should be 
mini-scale perturbations u" , v" , and p" comparable to co . 

The work of Molo-Christensen (1970, 1971) shows how these 
mini-scale perturbations influence dynamics of the large- 
scale motion. For the purpose of this document the complexity 
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of the double-prime scale is not needed, so the analysis 
retains only the single-prime scale plus the relaxation 
parameter to which is at a molecular scale or possibly at 
the scale of molecular clusters. 


Substituting Equations IV-9 into Equations IV-6, 7 gives 


*UL+*Ll + (U + u’) 9 < u ytn + (V + v’) 


Tt 

- _1 3(P + P’) , ^ 

p 5x A 2 


3x‘ 


(U + u') + 


3y' 


Ty~ 


(U + u’) 


+ A. 


(If) 


IV-10 


+ ( u * u') hv i-iD . ( v . v') ^ y > 


at 

.1 (P * P ') + A 

P ay 2 


2 2 

9 (V + V») + 3 (U + u') 
3x 3y^ 


- A. 


(If) 


IV-11 


If ■ 2B 3? - CU * 


“’>(&) - < v * V)(!f) • 


IV- 1 2 


It is clear (from Equation III-8) that the z-component of 
spin is 


l/3v 3u\ 

w " 7\77 ' 7y) 


(0 


1 

7 


(V ♦ V) - £ CU . 


u')J 


(0 


1 / 3V 3U\ 1 / 3v' _ 3u 1 \ 

7 \ 3x 3y / 2 \ 3x 3y / 

(ft + to ’ ) - (ft + <o' ) 


(0 


to 


to 


IV- 13 
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The term of 0 is identified as the vorticity of the mean 
flow, to' is identified as a perturbation in vorticity 
associated with the perturbations u' and v' , and 0 is 
identified as a mean relaxation parameter with an oscillatory 
modification given by to' . Even though the prime quantities 
are perturbations from the mean flow, they are of a completely 
different scale than to' , and the sum of 0 and to* may 
be taken as an equilibrium value to when considered 
relative to to . For an estimate of the scale of to' see 
Appendix A. 


Equation IV-12 then expands to 


8to _ 

7t " 


2B 3 (0 


(U + u') 


/ 30 
\ 3x 


+ 



30 3(o \ 

3"x TTx / 


- (V + v') 


/ 30_ 

\ 3y 


3to' 

w~ 


30 

3y 
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Substituting Equation IV-13 into Equation IV-10 and separating 
mean from perturbation terms gives 


3U 
3 1 


+ V 


+ T T 

9U + 

v 

3U ) 

( 3u 

+ U 

9x 

V 

3y J 

( 3 1 

3u ' 


3U 


3u ' 

3y 

+ V 1 

3y 

+ v* 

3y 


U 


3u 

3x 


+ u' 



u ' 


1 3JP _ 1 3p * 
p 3x p 3x 


3u 1 
3x 







3m ' 

^y~ " 3y " TT ) 


IV- 1 5 
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Doing the same thing for Equation IV-11 gives 

( 3V . IT 3V . v 3V) . Uv' . TI 9v ' . nf 9V . „ , 3v' . v 3v' 

U 37 V 3y iat - + u 35T + u 37 + u 37” V 3y” 


+ v * P + v* ^ 
3y 3y^ i 


P 3y 

2 


- I |El + A O 
p 3y 2)» 2 


1 3P 1 3p ' 
3"u' 


J OJ 


a"v + anj + a 2 v» 
ax" 3y 2 3x 2 

3ft aw* ) 


3y' 


I 9^ . 3( 

A 3|37 + 37” 3y ' TT 


IV- 16 


Now, paralleling the usual assumption that the mean flow by 
itself also satisfies the Navier-Stokes equations, the same 
assumption will be invoked here for modified Navier- 
Stokes equations with mean vorticity and relaxation terms. 
Thus, subtracting out the mean flow terms allows Equations 
IV- 1 5 and IV-16 to reduce to 

3t” U 37“ U 37 + u 37" V 3y” + v 3y v W~ 


3v 1 

at 



-I JEl + A I 

i a 2 u ' 

♦ M * 

A - l~' j 


p ax 2 
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ay 2 ) 

3 (ay ay ( 
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, av x . 
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,, av' 
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* U 8x 
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P IT 
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IV- 18 


where 


0) 


3v' 

3x 


3u 1 
3y 


IV- 19 
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The boundary layer assumptions are now invoked, so 

V ■ 0 ; - 0 ; If - o ; 

and since the perturbations are assumed small their products 
will be assumed negligible. It is emphasized that these 
several assumptions are not necessarily concurred with, and 
they may significantly influence the problem. However, they 
are made here in order to compare features of the structured 
fluid with the classical model used for hydrodynamic stability 
analyses. With the boundary layer assumptions, and dropping 
terms such as 


u' 


9u * 
9x 



etc . 
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Equations IV-17 and IV-18 reduce to 
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3 3 

Taking jy of Equation IV-22 and subtracting from it of 

Equation IV-23 eliminates the pressure terms: 
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It is now assumed that the disturbed motion is oscillatory* 
with the stream function given by 

iHx,y,t) = <j>(y)e i(ax “ ?t) IV-25 

where 

2 7T 

a = — = dimensional wave number 

£ = + if; = complex frequency IV- 26 

c = —■ = C r + ic^ = complex phase velocity 


F 

See Kistler 


(1971) 


? 
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and 


u' = = 4' (y)e 1 (ax Ct) IV- 2 

v' = -||- = -ia<Ky)e l(ax-?t) IV-2 

Note that primes on u and V indicate perturbation 
quantities, whereas primes on <}> indicate differentiation 
with respect to y . 

Then the necessary derivatives are 
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Inserting the proper derivatives into Equation IV-24 gives 
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This reduces to 
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To retrieve the classical Orr- Sommerfeld stability equation 
it is only necessary to let A 3 = 0 : 

(U - c)(<J>" - a 2 <J>) - <J>U" = -iH. [<)>»•• - 2a 2 <})" + ct 4 4>] . 

IV- 49 


Dividing all velocities by the maximum velocity of the 

laminar flow (or for boundary layers, use the free-stream 
velocity U^) and dividing all lengths by an appropriate 
reference length such as boundary layer thickness, the 
above equation is nondimensionalized to give 

(U - c)(<(>" - a 2 <j>) - cj>U M = [((>"" - 2ot 2 <t>" + a 4 <j>] IV-50 

which is the usual form of the Orr-Sommerfeld equation 
(Schlicting, 1955) with R as the Reynolds number. 


The two differences between Equation IV-49 and the more 
general equation for a structured continua (Equation IV-48) 
are the modified coefficient of the usual viscous part and 
then the additional terms in w * . Now the coefficient A 3 
may or may not be small depending on what it is assumed to 
represent (Appendix A). Supposing A is extremely small, 
as it would be if spin is taken as the spin of individual 

molecules, then the modification to the usual viscous solu- 
tion is negligible. However, the additional term 


A 3 ( 9 2 u)' 3 2 co 

ia ' 8x 2 8y 2 
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could have a strong influence on system dynamics even if 
is very small. This term may be written as 
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As with the classical theory, physical meaning is attached 
only to the real part, which is 
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From Equation IV-53 it appears that the physically significant 
part of the additional term is oscillatory and always damped. 
However, w' also will be oscillatory, and of the form 

, i (ax-qt) 

= y(y)e 
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Thus 
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and Equation IV-53 becomes 


! -a 2 y + Y"( e l(ax-nt) e 1 sin (ax - ( t) , 

/ 1 r 
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or 
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Equation IV-61 can be expanded and the real part is 


RE = 


a 


Y -aye 


2 _ i yw* 


cos (ax - n r t) sin (ax - S r t) . 
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Now it is clear that certain multiples of the spin frequencies 
n i may enforce the prime level perturbations and cause 
oscillations in the system to amplify. It should be noted 
that since n. is associated with molecular collisions it 
should be expected to be much greater than ^ and therefore 
able to drive this term to significant influence even though 
A 3 and y are small. 

While it is hard to believe that molecular level oscillations 
could feed much energy into macroscopic perturbations at the 
prime level, the work of Mollo-Christensen has clearly shown 
that there is actually an infinity of scales in between, and 
each need act only on its neighbor to produce a cascade type 
process in reverse. All theory aside, Mollo-Christensen 
also notes physical proof of the strong interaction between 
disparate scales of motion in nature. As one example fur- 
nishing dramatic proof, he was able to completely change the 
dynamics of the gravity waves on the ocean surface by covering 
it with thin surface slick which eliminated the tiny capillary 


waves . 
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V. Conclusions 

1. It is concluded from Section II that the constitutive 
equations for a structured continua can be obtained by 
proper interpretation of the velocity to which the stress 
tensor is proportional. This technique quickly yields 
the results of Dahler (1959) , Dahler and Scriven (1963) , 
and Condiff and Dahler (1964) , and does not have to 
resort to kinetic theory. 

2. It is concluded from Section III that the spin viscosity 
term which couples the system equations for balance of 
linear momentum and balance of intrinsic angular momentum 
is important since the coupling transforms the system 
from parabolic to parabolic plus hyperbolic, and allows 
the system to support vector wave propagation in the 
interior of the fluid in addition to the usual diffusion 
phenomena . 

/ 

/ 

3. From the development of the vector wave equation in 
Section III, it is concluded further that the system can 
be mechanically polarized due to flow. 

4. It also is concluded that classical sound waves can be 
coupled to the transverse waves. Thus, sound waves can 
excite transverse oscillations in the flow, and in turn 
the propagation of high frequency sound can be attenuated 
not only by classical dissipation but also by dissipating 
their energy into the transverse waves which they excite. 
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5. Since the constitutive and kinematic equations discussed 
in Sections II and III are representative of fluids such 
as air, it is concluded that air can support transverse 
waves . 

6. Since longitudinal and transverse waves in a physical 
media are usually called respectively compressional and 
shear waves, and since the ability to support both com- 
pressional and shear waves is representative of what is 
usually called visco-elastic, it is concluded that this 
work reasonably establishes a visco-eleastic character 
for air. 

7. From the stability analysis of Section IV it is concluded 
that the basic stability equation for the system has an 
additional term multiplied by a positive exponential term 
in time. Thus, a part of the viscous solution is 
inherently divergent and the system will be unstable 
unless it possesses sufficient damping in the other 
(i.e. the usual) part of the viscous solution. 

8. Since the unstable part has a very slow growth rate 

(i.e. 8 = small), it appears that in digital computations 

the methods used to suppress numerical instability would 
normally completely overshadow this physical instability. 
This would tend to imply that the instability due to 
skew-symmetric stress is a negligibly weak effect. 

Obviously, if the flow is tripped or has large free- 
stream disturbances, then the effects described herein 
are not going to influence the transition location much. 
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However, in cases of natural transition for which the 
free stream disturbances are very small and the viscosity 
is also small, then it is suggested that skew stress 
effects can be very significant in adjusting the actual 
transition point. 
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Appendix A: Order of Magnitude Considerations 


From Equation III-2, the case for no stress couples is 


A 

dw 


XT’ = |t (Vxu - 2w 0 ) » 
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which can be written as 


dw 

o 

cTT" 


■ ki- 


Vxu - to 


.) ■ 


1 ^ 

- (u 
T v e 


ft) ) 

o J 
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Note that this is like a simple reaction rate equation where 
/v X ^ 

w = Vxu is considered to be an equilibrium or relaxed 

G A T 

level of spin, and t = is the relaxation time. 


Since the relaxation is by molecular collisions, it is 
reasonable to conclude that the order of magnitude of the 
relaxation time is given by 


T 


Pi 

w 


0(collision time) 



0 / mean free path _ \ 
\avg. molec. speed / 


Now the ordinary coefficient of shear viscosity may be taken 
as 


P 


0.499 pvL , 
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Thus 


* - ft - HI) = °(*£) 

which implies that 

6 ■ °(H • 

In Section II, it was noted that in the theory of structured 

AA 

continua the I is a constant dyadic representing average 
moment of inertia of the rigid molecules in the flow. Thus 
Equation A- 7 be written in terms of the average radius 
of gyration of the rigid molecules, and this interpretation 

/\A 

of I would yield 


3 
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which is a very small quantity since mean free path is so 
much greater than the dimensions of a molecule. The coeffi- 
cient of spin viscosity, representing resistance to rotation 
of individual molecules, is therefore recognized to be quite 
small. Whether or not this small quantity and the effects 
it represents can be omitted is still open to question and 
must await extensive carefully analysis of the coupled model. 



There is second possible interpretation of I that must be 

AA 

considered. In Equation II-3, the parameter I actually 

is a sort of "moment of inertia of the collision", and for 

Equation A-7 this was interpreted to be the average moment 

of inertia of the molecules. Dr. R. C. Ried has noted 
, . /v\ 

another interpretation is that I could be an inertia 

tensor associated with molecular correlations or a kinetic 

theory formulation based on a doublet distributed function. 

In this case the average instantaneous radius of gyration would 

be approximately the mean free path. Hence k/L » 1 , and 

3 2 = 0 (y) , A-9 

which means that the spin coupling would be very significant 
in the problem. 

It is possible that both interpretations yield valid and 
significant terms, and that in Equation I I- 3 the intrinsic 
angular momentum must account for both the angular momentum 
change of a pair during collision plus the change of angular 
momentum of individual molecules about their own axis of 
rotation. 


The relaxation time associated with the first interpretation 
of I can be estimated as follows. For air at NTP, assume 


n 


a 


number density ~ 3xl0 19 
kinetic cross section x 
avg. molecular velocity 


molecules/cm' 


3*10 15 cm. 


5xl0 4 cm/sec, 
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v 



Then 


L = = mean free path 

~ 10 -5 cm A- 11 

and 

t = o/-\ = 0 /— — j) - 0(2xl0' 9 ) sec. A- 1 2 

W \5xl0 4 / 

This of course is a very short relaxation time. However, it 
is noted that propagation of ultrasonic waves in a gas is 
known to produce stress or temperature changes with a time 

_ 7 

scale on the order of at least 10 sec. For example, at a 
frequency of 30 Me, the period is 

t = i = = 3xl0 -8 sec. A- 13 

3x10' cy/sec 

Thus, high frequency noise radiating from a wind-tunnel 
boundary layer conceivably could excite such high frequencies. 
It is noted further that the flow under consideration will 
generate such high frequency waves as its own boundary con- 
dition, because at the wall the collision frequency is on 
the same order of magnitude as in the flow, and each colli- 
sion at the wall must be considered as collision with a 
rough surface and the generator of spin for the colliding 
molecules. It is postulated then that this boundary condition 
amounts to the excitation of the system with random disturb- 
ances and at frequencies of all magnitudes up to and including 
the maximum collision frequency. 



Section III showed that if the system can support the transfer 
of intrinsic angular momentum (via either spin of individual 
molecules, or spin of a colliding pair) then the rotation 
can propagate as a wave. It is thus postulated further 
that such spin originated at the boundary (or elsewhere) 
propagates into the flow and "pre-excites" local molecules 
so that they become more and more susceptible to clustering 
and random cluster rotation which could then be the physical 
mechanism of transition. 

AA 

It should be noted that if the I is interpreted as a 

moment of inertia of the binary collisions, then considerable 

/w 

care must be used in developing the analysis since I should 
be a function of time and other parameters. 
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Appendix B; Bulk Coefficient of Viscosity 

In recent years there has developed a controversy regarding 
the bulk viscosity of fluids. In some ways this has been 
healthy, because previously there was a widespread practice 
in classical fluid mechanics of accepting Stoke's relation 
as established truth. The matter is not a trivial academic 
argument since magnitude and functional dependence of any of 
the viscosity coefficients (y, X, or k) could have significant 
influence on stability of the fluid system (Kistler; 1971) . 

The manner in which these coefficients are determined experi- 
mentally (both as to accuracy and interpretation of the data) 
needs to be carefully re-examined with more attention on 
features important to a stability problem and not just with 
concern for thermodynamic considerations and computation of 
mean flows. 

For Newtonian fluids it generally is assumed that the elements 

M 

of D in Equation 11-14 are constants or functions of thermo- 
dynamic state, but that they do not depend on stress, rate of 
deformation, frequency, etc. However, in the field of ultra- 
sonics where dilitation effects are important, it now is 
customary to treat one or more of the viscosity coefficients 
as complex and frequency dependent as a means of adjusting 
hydrodynamic theory to fit data exhibiting relaxation 
effects (Tisza; 1952) . 

Truesde.il (1952 , 1953) questions the physical meaning of com- 
plex viscosity coefficients and of theory to fit certain 
data. Furthermore, he has shown that hydrodynamic theory 
cannot account for relaxation effects in general by a simple 
adjustment of bulk viscosity, because when it is adjusted 
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to make the theory fit absorption data in the linear range 
then the absorption is too great at higher frequencies. 

According to Truesdell, relaxation theory is only an ad hoc 
theory created to explain ultrasonic absorption and disper- 
sion. He holds that it is isolated from mechanics in general, 
and gives no idea how the fluid will behave for conditions 
other than infinitesimal oscillations. Truesdell claims the 
real implication of the theory is that fluids might have 
some purely mechanical properties not included in the classical 
hydrodynamical model given by the Navier-Stokes equation. 

This supports the intuition of Kistler (1969, 1971), and is 
in agreement with the equations of Dahler (1953), Grad (1952, 
1953) , and others . 

Now if such properties do exist, then they should manifest 
themselves in a variety of mechanical situations and not just 
in ultrasonics. Truesdell attributes this point of view to 
Gemant (1935) , who suggested that fluids exhibiting relaxation 
effects are in fact visco-elastic substances. Truesdell also 
indicates that formulae similar to the proposal by Tisza do 
result from the visco-elastic theories of Gemant and others. 
Furthermore, he has shown that when such formula are applied 
to plane infinitesimal waves, they yield an equation identical 
in form with the wave equation of classical elasticity theory 
except that the elastic moduli are replaced by certain func- 
tions of ico. 
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For the hydrodynamic stability problem it is likely that 
accuracy requirements on the coefficients are more stringent 
than those suitable for mean flow calculations, and any 
frequency dependence of the coefficients as already noted 
could be tremendously important. It is reasoned that experi- 
ments should be conducted up to very high ultrasonic fre- 
quencies because relaxation effects and mechanical polarization 
would be greatest under these conditions. It has already 
been suggested (Appendix A) that flow along the boundary pro- 
duces extremely high frequency waves (approaching the colli- 
sion frequency with the wall) , because each collision with 
the wall produces a translational and rotational perturbation 
to the system. 

Rosenhead (1954) notes that the significance of X does not 
enter into the usual phenomena of hydrodynamics. This is 
because it drops out of Equation 1 1 1 - 1 when the flow is 
incompressible. Furthermore, the acceptance of Stoke 's 
relation of course means that the bulk viscosity k is 
assumed to be zero, and it too does not enter the usual 
hydrodynamics problem. 

Regarding the coefficient X , when dilitation is signifi- 
cant (and even in cases where very small dilitation is 
responsible for phenomena of interest) , X cannot be neg- 
lected. This is well known to be the situation for absorption 
of sound waves in liquids and gases, and based on the results 
of Sections II and III of this document it is felt also to be 
the case for stability of laminar flows nearing transition. 

The coefficient X is usually taken as a constant connected 
with divergence of the velocity vector, but Andrade (1954) 
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points out that it is not a physical constant in the ordinary 
sense because no experiment can be devised which will measure 
it directly. In experiments where volume is compressed 
without change in shape, the coefficient involved actually 
is the bulk viscosity, k = (2y + 3A) . 

Rosenhead comments further that both y and A are intro- 
duced as quite independent phenomenological coefficients. 

It is only the assumption of k = 0 that leads to Stoke’s 

7 

relation A = y , which is equivalent to an assumption 
that the molecules have no internal structure. This is not 
generally valid, and actual fluids have values of A/y that 
may be positive rather than negative, with the magnitude of 
this ratio at times being as great as 200. 

It does seem to be reasonably well established now that y 
and A indeed are independent, but there still is considerable 
uncertainty about what should be used for bulk viscosity. 

The problem lies in the difficulty of devising experiments 
that clearly separate and distinguish influences of the 
various dissipative mechanisms. For example, Truesdell (1954) 
indicates that despite crude agreement between values of bulk 
viscosity calculated from absorption and from streaming, 
neither absorption nor streaming measurements, (in the linear 
range) can distinguish between the effects of various possible 
dissipative mechanisms such as bulk viscosity and non-linear 
shear viscosity. 
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According to Andrade, the streaming and sound absorption 
are interrelated rather than independent, and streaming at 
sonic frequencies seems to be a boundary- layer circulation 
due to sound absorption and explainable without resort to 
a second coefficient of viscosity. 

Clearly the whole problem of hydrodynamic stability is 
opening to a new line of thinking, and much work needs to 
be done in both the theoritical and the experimental areas 
to evaluate effects previously believed to be unimportant. 
Without significant experimental activities there is going 
to remain an inadequacy in the phenomenological data going 
into the mathematical models, and there will not be any 
suitable experimental results with which the analytical 
results can be compared. Without a rapid improvement in 
experimental activities there soon will be a gap very 
analogous to the gap between the theoritical work of Orr, 
Sommerfeld, Tollmien and Schicting, and the much later 
experimental confirmation by Schubeauer and Skramstad. 



